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1.0 INTRODUCTION 


The work which will be presented in this report was 
motivated by our observations that, if only the effects of 
dipole-dipole interactions are considered, theoretical calcu- 
lations of halfwidths for self-broadened H,0 transitions 
exhibit some rather serious discrepancies with recently re- 
ported tunable laser measurements .2222324/ These discrepan~ 
cies occur mainly at high and intermediate J values, and ob- 
tain using either the well-known Anderson~Tsao-Curnutte (ATC) 
theory228/ of pressure broadening, or a more recent method 


(QFT) developed by the present author. 228/ 


In Table I we present a comparison of theory with 
experiment for a number of high and intermediate J transi- 
tions. The results in Table I indicate that the theoretical 
widths are too small, with percentage errors sometimes ex- 
ceeding 50%. In most cases the results are only weakly de- 
pendent on which theory is used. However, all theoretical 
numbers are based on a calculation which includes only dipole-~ 
dipole interactions as the dominant scattering mechanism. 
These results suggest neglect of some important scattering 
mechanism at high J. Some possibilities are exchange scat- 
tering, dipole-quadrupole, or dipole-induced-dipole 
interactions. 


Our initial efforts to remove the above discrepan- 
cies have concentrated on the contribution of dipole-quadru- 
pole collisions between the H,0 molecules. A future journal 
article is planned in which detailed comparison of the ex- 
panded theory with experiment will be presented. In this re- 
port we will outline the computational methods necessary to 
obtain H,0 quadrupole moment matrix elements for the pressure- 
broadening calculations. This computation is far from trivial 
since H,0 is an asymmetric top molecule whose permanent quad- 


rupole moment cannot be represented by a single scalar quan- 


Comparison of Theoretical and Experimental H,0 
Self-Broadened Halfwidths 
(All halfwidths in cm '/atm; 

pure rotational transition; 


rot = 
+t measurement near 382°K, all others near 297°K) 


Theory 
Transition ATC QFT Experiment and_Reference 
8,1,8 + 9,4,5 (rot) “409 .402 -559, Ref. «u 
8,3,6 + 9,6,3 (rot) ~343 377 -611, Ref. 3 
11,2510 + 22,4,9 {rot} <261 «257 -319, Ref. 3 
L2¢2@yh]. + 1353522 Vo -160 .2140 -228, Ref. 2 
-170, Ref. 4 
15§,1,15 + 16,0,16 Vo -098 .083 Lit, Rete. 1, 2 | 
-166, Ref. 4 
38,2,14 + 16,1,15 Vo -084 .085 -203, Ref. 2 
17,0,17 + 18,3,16 (rot) .085 .078 -125t, Ref. 3 
19,0,19 20,3,18 (rot) .063 .063 -O89t, Ref. 3 


tity. The theoretical framework for the calculations is de- 
veloped in Section 2.0 and Appendix A. Numerical results for 


H,O quadrupole moment matrix elements are presented in Ap- 


2 


pendix B for some low J quadrupole transitions. 


2.0 PERMANENT QUADRUPOLE MOMENT MATRIX ELEMENTS FOR 
ASYMMETRIC ROTOR MOLECULES 


The calculation of pressure broadening due to quad- 
rupole interactions can be expressed! completely in terms of 
an m-independent, reduced matrix element, <J{{Qt|{a'>, of the — 
traceless quadrupole moment tensor Cc. the latter defined as 


t 


1 
} e; (Pr, ¥, - 37,’ rt. (1) 


g - 5 tr(g) f, co) 


where the sum is over all charged particles in the molecuic . 
and where f is the unit tensor, such that f.4 = Kef = K fo 
an arbitrary vector A. In Ea« ‘(€2), $ is the quadrupole me ee. 


tensor 
> > . 
§ = } e, T; Py> (3) 
and we note 
tr(§*) = Tr(Q) - 3 Tr(G) Tr(f) = 0, (4) 


since the trace of the unit tensor equals 3. 


The reduced matrix element <3{ {qt} fu'> is defined 
using the Wigner-Eckhart theorem2212 for traceless symmetric ta 
tensors, via the equation 


<J m(O"|u! m'> 


= <g}Jat|lu'> on a p (J'2-m'=-M|g'2 J-m) #,, (5) 
ee > 
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taken to be the direction of quantization for the magnetic 
quantum numbers m, m'. For non-linear molecules |J> actually 
stands for all quantum numbers (except m) which are necessary 
to specify the state, i.e. [J> = |[Jtee>. 


The reduced matrix element is then obtained by calcu- 
lating any component of <J Tt m{O*\a' t' m'>, e.g. the Z Z 
component, which is diagonal in m, m'. Thus, we have 


t tot 
<J tilet}fa' t'> = Spe eas : (6) 
(J'2-m O|J'2 J-m) 

The next step is to transform to a principal axis 
system (x,y,Z) which diagonalizes ét and 6. In the case of 
HO, it is clear from symmetry that the principal axes for 
the quadrupole moment and inertial tensors are identical 
except possibly for the choice of origin. However, in most 
pressure-broadening theories this choice is fixed to be the 
center of mass. This results because one typically writes 
the various multipole interactions as V(R;-R,) where Rs R, 
specify the centers of mass of the two interacting molecules. 
This point is somewhat subtle and is discussed more fully in 
Ref. (17). The transformation to the principal axis (x,y,z) 


system can be written 


t 1 
<J t m|@75[0' <' m> = ee : a m| C4. - x] |J' t' m>, 
one > (7) 


= , fe. Ps, var inci i 7 
where Qais } e; Tig Pia e the principal axis components of 


6, and where 4 = Zea are the direction cosine operators 


connecting the two coordinate systems. 


The matrix elements of (O54 - #1 can be eva! ated 
by transforming to a representation |J K m> of symmetri-_ top 
eigenstates. Some further details of this evaluation are 
given in Appendix A. 


The result for the reduced matrix element 
<J ti lat} fa" t'> is then found to be 


sf xlie lata =1Gir 
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*4(Q a le & Quy 1g" «te 


- js | <3 4d K> A 
ZZ KK! 


ceent 
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3 - ' ' ' 
-Vr (Gs hy) bs <0 2/0 WB srg <o* Ka 
where <J' K'|J' t'> are the eigenvector components of state 
|J' t'> in the symmetric top representation. The A and B 
matrices in Eq. (7) are given explicitly by 


A = js (o* 2 et elgt 2 a KIS 
JK3J'K' ~I+] K,K"? 


1 f2g'+1 


Paxso'k! ~ Foy 2041 


{(o' 2 Kt 2]gt 2 TCK'#2)) by vagy 


# to! 2 Kt = 2]d* 2 GOC'H2)) 6 Le 


K,Kt<2 


and it may be verified that A, B satisfy the sum rules 


A B 


*IK3I°K! *IK;g'K! = Ts (10a) 


} } 
o'K! g'K' 


pa Axxsgtxt Boysgtkt = 2% (10b) 
Several points should be noted before going on. 
First, from Eq. (8) we note that the reduced matrix element 


for an asymmetric rotor depends on two independent scalar 
parameters 


Q, = (Qi - Fy e (11a) 


Qo = (,. - & + (11b) 


However, this parameterization can be expressed in other use- 
ful fashions. From Eq. (2) we note that 


ac At 


Qi 7 = z Qyy?? (12a) 


Q, = (Qe, = (12b) 
where the atts are the principal axis components of the trace- 
less tensor a. These three components are not independent 
: t t ts 
since es + Qyy + Oe = 0. Hence Qy an can be parameter- 
ized in terms of any two components of Q°. 
Secondly, it should be pointed out in Eq. (8) that 
the labels (x,y,z) refer to an arbitrary symmetric top repre- 
sentation (arbitrary, except that (x,y,z) must be principal 
axes). In particular, z need not correspond to the (two-fold) 
symmetry axis of the H,0 molecule. In fact the numerical re- 
sults which we present in Appendix B are based on a IR repre- 
sentation, where (z,x,y) are associated with the lowest-order 
constants, A > B > C, in the rotational-vibrational Hamiltonian 
via 


g 2 2 2 
H = B ve +c Jy + A Jo (13) 


For H,0 the intermediate constant B represents the inverse 
moment of inertia about the symmetry axis. 


Of course for a symmetric top molecule, it is most 
natural to choose z, the symmetry axis of the symmetric top 
representation, to coincide with the symmetry axis of the 
molecule. In this case ‘oy - Qyy? = 0, the second term in 
Eq. (8) vanishes, and the equation collapses to give 


<3 t{I1Q*}]u* ct> 


Jit 


= <3 K[[Qt{[gt Kt> = (41 “1y) 


Qy AgKsatKt? 


where Q, = (Q,. - Qe) corresponds to the Buckingham’ , {| 
Birnbaum— definition of the scalar quadrupole moment for 
symmetric tops. This definition is one-half the value used, : i 
e.g.- by Tsao and Curnutte®/ and Benedict and Kaplan23/ in 4 


discussing pressure-broadening. 


For the purposes of numerical computation it is 
advantageous to transform Eq. (8) further using tc Wang 
transformationdt2lS/ | which corresponds to formin, i 
combinations of the symmetric top wave functions accord..._ 


[E,(0 K)> = Lf Ja K> + [u-K>} K even > 0, (fa) 

"2 | 
= |J O> (K = 0) 

[EG K)> = Sede [oe to # fae) K even > 0, (15h) 
x? 

0,0 Ki = Zod lg ee + Joe) K odd > 0, (15c) 

+ 

vo 

[o_(I K)> = a lJ K> + |g-K>} K odd > 0. (15a) 


Let us write the above Wang (symmetric top) states as |J K v> 
where v = 1, 2, 3, 4 = Ey, E_, O,, 
metry. The exact eigenstates [J 1t> of the asymmetric rotor 


O_ denotes the Wang sym- 


can also be classified according to their Wang symmetry, i.e. 
JJ t> = [J t v>. In this notation Eq. (8) can be written 


> 


‘ 
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3 
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q 
< J <8 « vid K ve B K* v*l|g" ef yte}, (16) 


IKv;giKtyt <I! 
Here, e.g., 


A 


ekegiery = <OLAE WlaloAst a> 
Sieg 


= % (~ <I K] + <J-K]) A (- [ot K'> + [ut-K'>) 


ick 
= 5 {Asya * Agexsat-x) 8x,x" 

S = afk 

= Asyigrx 8x,xr for AJ = J - Jt = 0, #2 (17a) 
< § for Ad = #1. (17b) 


The results for all other A and B matrix elements can be ob- 
tained in a similar fashion, and it may readily be verified 
that the non-vanishing matrix elements of both A and B obey 
the selection rules 


E, © £5 EL * Bis 0, © Dyy 0, © 0, for ke = 0, £2, 
E, + E_, 0, + 0_ for AJ = #1. 


That there are no allowed matrix elements between even and 
odd states follows from Eqs. (9) since AK = K - K' is always 


even. 
We now define 
Agaysa'tty! i 
Yo <3 t v[J K v> As serge SOP KI vifgl tt vio, (18a) 
KK f > 
9 


liberate eer 


Borvst't'y! 


<3* Kt gthg® ¢* gts, 18b) 


une <J t v[J K v> BoxysatKty! 


These matrices satisfy sum rules similar to Eqs. (10), in 
particular 
aes (19a) 


2 = 2 
AScusd't'y! E I BItusd' tty! 


J't'y!' J'tty! 


Jtztyt Agcysa'tty! Burysatrtyt = 9 (19b) 


In terms of these matrices, Fq. (16) assumes its firal forr 


<I t v{fat{{gt tt vt> = aie 


t i AE dat 
; {(Q7,, — FM 7 F Qyy? Agrysattty! 
S: t t 
~ % Qe > Qyy? Brrvsatatytts le 


Programs to compute the A, B matrices for H,0 ana 
other asymmetric rotor molecules have been written and are 
presently operational. These programs operate as subroutines 
which are attached to programs, developed at AFGL, which gen- 
erate the eigenvalues and eigenvectors of the asymmetric 
rotor Hamiltonian. These latter programs also include vibra- 
tional-rotational coupling via the Watson Hamiltonian method. 
The programs are quite efficient, requiring approximately 60 
seconds of CP time to compute the strongest (AK, = 0, +2) 
quadrupole transitions for J < 22. A program listing of our 
quadrupole moment subroutines is furnished in Appendix C of 
the present report. 


Numerical results for the A, B matrix elements in 
Eq. (20) are presented in Appendix B for some low J quadrupole 
transitions. The labeling of the K,> Ko states in the tabu- 
lation of Appendix B is consistent with a IR symmetric top 
representation for H,0. We have verified that our methods do 


give representation-independent results, in the sense that 


the same matrix element <J tllot{ fa t'> always corresponds 
to the same energies Esq? Exret and energy difference Est - 
Esiqte However, if one switches representations, e.g. from 
IR to IIR, the K> K labeling of the states also switches, 


as discussed in Allen and Cross. 28/ 


Finally, to compute <J t{|Qt||ug' t'>, one needs to 
know values of Cae Vy? 2 (which, again, are not all inde- 
pendent since Tr(Q-) = 0). For H,0 there appears to be no 
experimental determination of these quantities, however, 
Stogryn and Stogryntl/ have listed fairly consistent theo- 
retical values from quantum mechanical calculations. Their 
values are also quoted relative to the center of mass of the 
molecule. From Table 3 of Ref. 17 we deduce the following 


average values: 


Ge = .12 x 10°?® esu-cm? 
a “26 2 

Qyy = -.72 x 10 esu-cm 

uc. = .60 x 10 7° esu-cm? 


In order to make a simple comparison with the No 
quadrupole moment, we note from Eqs. (19), (20) that 


2 
I [<a tilat}}at rt>| 
J't* 


t a Lats” 4 Seat < 
zz ~ 7% ~ F yy? + FT Que - Qyy? 


(Q 


= 1.34 x 10°52 (from values listed above). 


The corresponding result for No is 


2 
Y |<affatjat>|? = (1.52 x 10728 esu-cm?) 
J' 


— 3 er --. 


» 2.38 «0 **, 


- i asi as ee si lait io ais tee ee EST HY ae areca aires Mis ee 


en ne a eet a nr Semen emt nema ee ~ See noi MOI — 


where the numerical value has been taken from Table 1 of 
Ref. 17. 


Since the strength of dipole-quadrupole scatter_ng 
in pressure broadening theory is directly proportional to the 
quadrupole moment squared, we infer from the numbers given 
above that the strength of this type of scattering for H,0- 
H,0 should be roughly 60% of the corresponding strength for 
H 
tions that the H,0-N, dipole-quadrupole strength is typically 
20% of the H,0-H,0 dipole-dipole strength. On this simple 
basis we expect the H,0-H,0 dipole-quadrupole collisions to 


29-No- In turn, one knows from pressure broadening calcula- 


contribute something like a 10% correction to the 
broadened halfwidths. 


The effect of these collisions will sometimes be 
enhanced when resonant energy denominators occur in the 
theory, however, this enhancement cannot be estimated witho t 
performing detailed calculations. Preliminary results fra | 
our pressure broadening programs indicate that the dipole- 
quadrupole collisions for H,0-H,0 are negligible at low J, 
while at high J (J > 12) they can contribute 10-15% of the 
total linewidth. This suggests that other mechanisms such 


as exchange scattering need to be considered. 
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APPENDIX A 


Shakes 
Matrix Elements of cee. *] 


We wish to evaluate Eq. (7) by introducing a com- 


plete set 


of symmetric top states. 


y jd Km <d Km{ = 1 
Km 


Since J, m are also good quantum 


numbers for the asymmetric rotor, we obtain 


<J tT m[Qols' 


ee lee 


O=X,y 5Z 


The symmetric 


Vaqmies¥2o) = 
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The 
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which are independent of 6. 
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(o. m 73) ae 
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<J' K' mijg' +c! m>. (Al) 

top wave functions have the form®’ 
imd ,ikKy 

e e f sxm‘®? » (A2) 

are the three Euler angles. 


direction cosines (Oo .9 bry? > can be ex- 


cosé 


sin@sinyp 


sin6cosp : (A3) 


It follows that 


e Y,968.¥) 5 (Au) 


ee 


RR 


‘ un 
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~ [FE (%p9 (OW) + Yo (0,H)) 5 (AS) 


At y 
Cay = 3) = = FV Ee Log (Ov) 


+ [Fe (X59 (850) + YX, (09)? (A6) 


where the Y's are spherical harmonics with the Condon and 
Shortley choice of phases. 


The matrix element which Pi needed in Eq. (Al) are 
readily obtained using the formula”! 


<J K mY.) (8,0) |J' K' m'> 


= (- ay /2eeh jae (IH AKT [TY 5 TOK) 6 jays 


- (J' j m' Olu' j J m') Sum (A7) 
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MATRIX ELEMENTS OF A AND B (EQS. 18-20) FOR HO INA 
IR REPRESENTATION FOR SOME LOW J QUADRUPOLE TRANSITIONS. 
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APPENDIX B 


Matrix Elements of A and B (Eqs. 18-20) for H,0 ina 
IR Representation for Some Low J Quadrupole Transitions. 
(A and B elements are dimensionless and satisfy the 
sum rules in Eq. (19); energies are in units cm '.) 


J, Ki, K, J's Ki, K Es E,; 7 Es, <JJAlu'> <dJ|Blu'> 
o, 6, 6 2. 8, 2 00.00 - 70.09 - 99040 .13822 
6. te-8 2. 2, 6 00.00 -136.16 -.13822 - 99040 
a i 2, Gy 2 23.79 00.00 -.63246 .00000 
t, G5 3 ae ee | 23.79 <1i1.21 .00000 -.81650 
t, 6, a, 6, 3 23.79 -112.97 -74391 .16089 
1, G2 Sy 25 2 23.79 -188.36 -.21586 55448 
5 a Se se 42.37 00.00 «31623 -.54772 
|e ee | 2° 24-2 42.37. = 37.13 70712 ~40825 
Ls bs 6 5, se 2 42.37 -130.99 .62988 -.11805 
¥, hy 0 8, 8, 0 42.37 -243.05 -.05704 . 72069 
ti ts 3 iy dy 2 37.14 00.00 31623 .54772 
ay tet ae ee | 37.14 - 58.04 TOTAL -.40825 
dy Ay a 3, Es 3 37.14  -105.14 - 63088 . 23203 
2, 3.3 ay 3, 1 37.14 -248.08 -. 04465 -69245 
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PROGRAM LISTING OF QUADRUPOLE MOMENT SUBROUTINES 
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